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We derive an exact solution of an explicitly time-dependent multichannel model of quantum mechanical 
nonadiabatic transitions. In the limit N 1, where N is the number of states, we find that the survival 
probability of the initially populated state remains finite despite an almost arbitrary choice of a large number of 
parameters. This observation proves that quantum mechanical nonadiabatic transitions among a large number 
of states can effectively keep memory about the initial state of the system. This property can lead to a strongly 
non-ergodic behavior even in the thermodynamic limit of some systems with a broad distribution of coupling 
constants and the lack of energy conservation. 



Quantum mechanical nonadiabatic transitions have been 
studied in molecular collisions for long time CHS. Relatively 
recently, multichannel nonadiabatic processes became in the 
focus of research in condensed matter physics due to the in- 
terest in control of quantum many-body states in electronics, 
magnetic systems, and Bose condensates EHT). Condensed 
matter physicists often achieve an insight about their systems 
by finding a proper model from a large variety of exactly solv- 
able systems with well defined thermodynamic limit. Unfor- 
tunately, most of such known models correspond to the sta- 
tionary Schrodinger equation. Exact solutions of explicitly 
time-dependent systems with a large number of states are rare. 

The multistate version of the Stueckelberg-Majorana- 
Landau-Zener (LZ) model is one of the most investi- 
gated explicitly time-dependent problems fTQlfHl . Its goal 
is to find transition probabilities among N discrete states in- 
duced during the time evolution from r = -oo tor = -f oo 
in systems described by the Schodinger equation with linearly 
time-dependent coefficients: 

i^ = (A + BT)i>, (l) 

where A and B are constant N x N matrices. Matrix B is 
diagonal. Eigenstates of B are called the diabatic states and 
off-diagonal elements of the matrix A are called the coupling 
constants. One has to find the scattering N x TV matrix S, 
which element S nn > is the amplitude of the diabatic state n' at 
r — » -f oo, given that at r — » — oo the system was at the state 
n. The related matrix P, P nn > — \S nn > r, is called the matrix 
of transition probabilities. A number of exact results in which 
the regime N ^> 1 can be effectively considered, such as the 
Demkov-Osherov [12 ] model, are currently known. 

Available exact results seem to indicate that transition prob- 
abilities in the model ([I]) behave somewhat "classically": All 
known exact solutions with a finite number of states can be 
either reduced to independent 2-state systems lfl3lL or can 
be completely understood in terms of individual pairwise LZ 
transitions that happen in a chronological order between each 
pair of diabatic states (TUJ [121 031. Even beyond Eq. ([I]), 
Ostrovsky studied a similar model to ([T]) with nonlinear level 
crossings and came to similar conclusions fT5ll . 

Numerical simulations show deviations from the semiclas- 
sical behavior in non-integrable cases but the available exact 




FIG. 1. Time-dependence of the eigenspectrum (adiabatic energies) 
of the Hamiltonian with elements Hoo — k 2 /r, H nn — f3 n r, 
H n o — Hon — g n , and zero otherwise (n = 1, . . . , N; N = 8). 



results and numerical simulations of the model ([T]) with large 
sizes of the phase space lfT6l suggest that dynamics of such 
driven multistate systems can be at least qualitatively and ap- 
proximately understood in terms of the successive two-state 
interactions that happen at moments when pairs of diabatic 
energies cross or pass close to each other. Pairwise transition 
amplitudes then can be used to estimate amplitudes of the full 
trajectories in the phase space. 

Based mostly on the available exact solutions of ([T]), it has 
been argued in various contexts 1 7 ] that if a level with an ex- 
tremal slope crosses diabatic energies of a large number N of 
initially unpopulated states then the probability Po^o to re- 
main in the initial state after all intersections decreases expo- 
nentially with N, i.e. Po^o ~ exp(-aN), where a is a coef- 
ficient, which characterizes a transition probability at a typical 
encountered avoided crossing point. Such a result would be 
perfectly justified in classical stochastic kinetics ifTTl . Quan- 
tum mechanical exact results only extended the validity of 
such semiclassical expectations. There can certainly be excep- 
tions that correspond to resonances that can encounter when 
a large number of parameters are fine-tuned but this does not 
change the general qualitative understanding. 

In this letter, we provide an explicit example, which un- 
ambiguously demonstrates that the above intuition about the 
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FIG. 2. (a) The integration contour A incloses all branch cuts 
(dashed lines) from a large distance. The integral over this contour 
can then be expressed as the sum of the integrals over contours ~f n , 
n — 0, . . . , N along the branch cuts, (b) By a change of variables, 
each integral over j n can be transformed into the integral over the 
contour C. 



large- N limit of quantum mechanical nonadiabatic transitions 
becomes incorrect if we consider the evolution beyond the lin- 
ear time dependence of diabatic energies. We will show that 
even if N is large and all parameters in the Hamiltonian are 
chosen almost arbitrarily, e.g. from some random distribution, 
the probability to remain at the initial level after all level in- 
tersections can remain 0(1), i.e. essentially finite. 

To achieve our goal, we will solve a model of (N + 1) 
interacting states with the evolution equation for amplitudes 



N 



i — a 
dr 



dr 



PnTb n + g n a, (2) 



where n = 1, ... TV. We will call the diabatic state with the 
amplitude a the 0-th level. Its diabatic energy (i.e the energy 
that it would have at zero coupling to other states) decays ac- 
cording to the Coulomb law, ~ 1/r, while diabatic energies of 
other states change linearly, as shown in Fig. [T] For this prop- 
erty, we will call the model ^ the LZC-model, after Landau, 
Zener and Coulomb. At r — )> 0, transitions among states do 
not happen because the 0-th state has infinite energy. We will 
assume that the 0-th level is populated at r — ^ 0, while other 
states have vanishing amplitudes. During the evolution, for 
some time, distances between diabatic energies become com- 
parable to coupling constants. This is the time when nonadi- 
abatic transitions mostly happen. At r — > +oo all diabatic 
energies become again well separated and transitions termi- 
nate. 

From the semiclassical picture, we expect that the behav- 
ior of the transition probabilities will be mainly influenced by 
the encountered crossing points of the 0-th level with levels 
n having f3 n > 0, as shown in Fig. [T] Drawing the intuition 
from [7], one may also expect that the probability to remain 
at the initial 0-th level after all intersections will be exponen- 
tially suppressed at large N. 

We will solve the LZC model by changing variables as a = 
rbo and t = r 2 /2, leading to 



N 



dt v 



We then introduce the anzatz 



71=1 



.dbn f3 n b n ^g n b . (3) 
dt 



»(*) = / 

J A 



due- lut B n (u), n = 0, 



(4) 



where A is a contour such that the integrand vanishes at this 
contour limits. Substituting ^ in ^ we obtain a 1st order dif- 
ferential equation for Bq(u), which is trivially solvable. Sub- 
stituting the result back to ^ we find 



bo(t) 



'L 



due- mt (-u) 



N 

n 



-u + p ri 



(5) 



bj(f) = -Qg. 



L 



du ■ 



N 

n 

n=l 



-u + p n ye 



Integrals in ([5]) have almost the same structure and the same 
initial conditions for amplitudes as in the model of Ostrovsky 
ifTSTl . Following his steps, we consider a contour A that in- 
closes branch cuts at u = ({/? n }> 0) fr° m a large distance and 
goes to infinities at u = —ioc±R, where R is a large real 
number. In this limit, we can disregard terms /3 n in compari- 
son with u, so that integrals in ^ simplify, e.g., 



bo(t) Q I e 

I A 



-iut 



(-«)• 



2 du. 



(6) 



In ([6b, the contour A can be transformed into the contour C in 
Fig7[2jb) by switching to the variable z = iut and shrinking 
the contour to run around the branch cut of z. We can then use 
the formula 



T(z) 



2ism(7rz) J c 



{-r) z - l e- T dr 



(7) 



to evaluate the integral. Returning to original variables: 
a{r) = rb (r 2 /2) and b n (r) = 6 n (r 2 /2), we find that, at 
r = 0+, all b n vanish, while a remains finite. If we set 



+3i7r/4-iar g r(l/2+ifc 2 /2)-ifc 2 ln(2)/2 



then we will arrange at r — » + that 



irk 2 



(8) 



a(r) r ^ 



~ r 



-ik z 



^•(r) r ^o-0(r), j^O, (9) 



which corresponds to the initially populated 0-th level asymp- 
totic. In order to find transition probabilities at r — > +oo 
limit, we continuously deform the contour A into a combina- 
tion of contours 7 n that inclose the branch cuts at u = f3 n as 
shown in Fig. [2ja). In the limit t —> +oo, only the vicinity 
of the branching points contribute essentially to each integral 
over 7 n . Hence one can change variables u — » u + /3 n , keep- 
ing the dependence on u only for terms that are singular near 
the origin of the 7 n -th cut. In all other factors, we can substi- 
tute u by its value at this point. Simple dimensionality argu- 
ments show that the integration for bj(t)t^.+oo in ^ over 7^ 
remains finite only if i = j . Hence a j-th integral in ([5]) over 
7j provides the asymptotic at t — » +00 for bj(t), i.e. 
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where we should assume that (— i) = e~ zn ^ 2 and —1 = e~ Z7r . 
Remaining integrals again can be evaluated with Eq. ([7]). In 
the supplementary file lfT8l we provide explicit expressions 
for asymptotic amplitudes at r — » +oo that we obtained 



from ( |1Q| ). Transition probabilities can be obtained by taking 
squares of the absolute values of transition amplitudes. Sur- 
prisingly, the final result appears very simple looking. In order 
to write it, it is convenient to introduce LZ-like probabilities: 



Pj 



j = l,...,7V. 



(11) 



The transition probabilities from the initially populated 0-th 
level to all possible states are then given by 



11 Pn 



„ //3m<0 

II Pn 



0^0 



1 + e 



-tt/c 2 



(12) 



II Pn (1 - ft) 



ft > : P ( 



0->j 



1 



t/c 2 



(13) 



II (l-Pi)c-* 



ft- < : P ->j 



1 + e 



-TT/C 2 



(14) 



/3n>0 



Here, for example, [| Pn means the product over all p n with 

n 

indexes n of states having slopes f3 n larger than zero. If there 
are no such states, this product is set to unity. Several qualita- 
tive observations follow straightforwardly: 

Corollary 1: The transition probability from the 0-th level 
to the level j with ft > depends only on coupling constants 
and slopes of the levels n with f3 n > ft . 

Corollary 2: The transition probability from the 0-th level 
to the level j with ft < depends only on coupling constants 
and slopes of the levels n with /3 n < ft . 

Corollary 1 would be expected from the semiclassical point 
of view. In the independent crossing approximation, only 
crossing points that precede the crossing of the level j would 
influence the population of the j-th level. This situation 
has previously encountered in the solution of the Demkov- 
Osherov model fV2\ . Parameter k characterizes the curvature 
of the 0-th level near the avoided crossings lfT8l . The limit of 



large values of k corresponds to locally linear dependence of 
0-th diabatic energy on time, so that the formulas ([T2|)-([T3]) 
naturally reproduce the results of the Demkov-Osherov for- 
mula 1 12 ] if we set k —> +oo. 

In contrast, Corollary 2 is very counterintuitive. Surpris- 
ingly, none of the encountered avoided crossing points influ- 
ences transition probabilities to the levels j with ft < 0, no 
matter how strong can couplings at crossing poins be. More- 
over, only couplings to the levels n with f3 n < ft , i.e. the lev- 
els that diverge faster than the level j from the 0-th level, con- 
tribute to the transition probability to the j-th diabatic state. 

In Fig. [3] we compare our theoretical predictions ( 12 )-( 14 ) 
with transition probabilities obtained in our numerical stud- 
ies of the evolution ^ in a three- state LZC model. We find 
an excellent agreement between numerical and theoretical re- 
sults. In particular, Fig. |3ja) shows that if ft > ft > 
then the transition probability to the level with the slope ft 
does not depend on the coupling #2 to the level with a lower 
slope, in agreement with Corollary 1 . In case when ft > ft 
and simultaneously ft < 0, we find that already the transition 
probability to the level with the slope ft does not depend on 
the coupling gi between the 0-th level and the level with the 
slope ft . 

Corollary 3: The probability to stay at the 0-th level at 
r — » +oo remains finite if either ft > or ft < for all 
j = 1, ... TV. Moreover, in the limit N — >• oo and at finite 
couplings, the probability to remain in the 0-th level depends 
only on the parameter k. 

The last observation is the main result of this letter. It is 
normally expected that if a 0-th level interacts with a large 
number of states with arbitrary coupling strengths, such an in- 
teraction would lead to a sort of equilibration. It has been 
often argued that the survival probability should be, actu- 
ally, vanishing exponentially with the number of levels N 
that interact with the highest slope level |7|. The 0-th level 
satisfies this condition if conditions of Corollary 3 are sat- 



isfied. Nevertheless, Eq. (12) predicts that if all ft > 
then P ^o > e'^ j{\ + e" 7 ^ 2 ), and if all ft < then 
Po^ >l/(l + e- fc2 ). 

In conclusion, we determined transition probabilities in a 
multichannel model of nonadiabatic transitions with a well 
defined transition probability matrix among all diabatic states. 
In our model, one of the diabatic energy levels changed non- 
linearly, which allowed us to investigate the effect of the cur- 
vature of this level on the transition probability matrix. We 
found that the behavior of the probability to remain at the ini- 
tially populated level, as well as the levels having a negative 
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FIG. 3. Probabilities of diabatic states at time t - 
LZC model for (a) ft > ft > 0; k 2 = 0.7, g 1 
ft 0.5; (b) ft > 0, ft < 0; k 2 0.5, <? 2 
ft = -1; (c) ft < 0, ft < 0; k 2 = 0.2, <? 2 = 0.3, ft - -0.5, 
ft = —1. Solid curves are theoretical predictions by Eqs. |T2|-|T4|. 
Discrete points represent the results of the numerical calculations. 
Insets on the right illustrate time-dependence of diabatic energies. In 
all numerical tests, the time evolves from r = 0.0001 to r = 800. 



slope of diabatic energy, can behave very counterintuitively. 
For example, some of the transition probabilities can be finite 
but independent of the couplings at all avoided crossing points 
that encounter during the system's evolution. 



We predict that systems with nonadiabatic transitions in the 
large N limit can show strongly non-classical behavior, e.g. 
with some states keeping a finite amplitude despite a large dis- 
tribution of coupling constants. In supplementary file |[T8lL we 
demonstrate numerically that this property appears also in a 
wide class of diabatic potentials, e.g. of the form ~ l/t r with 
some constant parameter r. The independent crossing approx- 
imation works for them only when Po^o remains close to 
unity, and the exponentially strong enhancement of the prob- 
ability to stay on the initial level is observed otherwise. 



The Coulomb potential frequently appears in studies of 
nonadiabatic transitions in molecular physics d [15] H9l . 
Moreover, if one of the levels is degenerate before a molecular 
collision or an application of a time-dependent field, this de- 
generacy experiences fan-like spitting, as in Rydberg's atoms 
due to the Stark effect |20|. Hence, the LZC-model repre- 
sents an important class of processes in atomic and molecular 
physics. 



Among other applications, we mention that, in the limit 
k — )> 0, the LZC-model corresponds to the previously unstud- 
ied case of an N-state bow-tie model lfT4l with the time evolu- 
tion from r = to r = +00. In the supplementary file lfT8ll . 
we show that the two- state version of the LZC model pro- 
vides, at a little cost for complexity, a better approximation 
of the transition probability at an avoided crossing than the 
Stueckelberg-Majorana-Landau-Zener formula. It also pro- 
vides a sample of a solvable model of an avoided diabatic 
crossing with a constant off-diagonal coupling. Finally, we 
provide an example in [18 ] of a possible realization of a 3- 
state LZC-model in an interacting qubit system. 
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SUPPLEMENTARY MATERIAL 

In support of the main text, section 1 of this supplemen- 
tary material provides the explicit form of the asymptotic 
values of amplitudes. Section 2 explains details of the nu- 
merical algorithm and provides additional numerical proof 
of the theoretical formulas in five-state LZC models. Sec- 
tion 3 discusses similar counterintuitive behavior in noninte- 
grable models with a singular potential of 0-th level. Section 4 
shows an example of the two-qubit Hamiltonian that realizes 
the three-state version of the LZC-model. Section 5 explores 
the 2- state version of the LZC model and its applications as 
an alternative to the LZ-formula and as a model for a diabatic 
avoided level crossing. 



I 

I. TRANSITION AMPLITUDES IN LZC-MODEL 

Asymptotically at r — >> 0, we choose amplitudes so that the 0-th level is initially populated. We fix constant phases so that 

a(r) T ^ ~ r~ lk \ yrU-O(r), j = l,...N. (15) 



The amplitudes at r — >> +oo are then given by 
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(20) 



II. NUMERICAL TESTS 

In order to check our theoretical predictions, we wrote a nu- 
merical program that simulated the quantum mechanical evo- 
lution that starts at a small initial time moment t = dt and 
proceeds in discrete time steps of sizes equal to dt. The evo- 
lution at each time step is described by the evolution operator 

U(t) = (l + iH(t)dt/2^j (l - iH(t)dt/2) * , (21) 

where H(t) is the time-dependent matrix Hamiltonian, and 
1 is the unit matrix. Operator U(t) is unitary. Up to 
terms o(dt 2 ), it coincides with the true evolution operator, 

l exp (ifi +dt H(t)dt), where f is the time-ordering oper- 
ator. Being equivalent to the true evolution operator in the 
dt — ^ limit, operator U (t) is relatively easy to calculate nu- 
merically, because the hardest step is the calculation of the 

inverse of the matrix ^1 — iH{t)dt/2 s j . By acting with U (t) 

on the state vector, the latter was updated up to large times at 
which no changes of the state probabilities was observed up 
to two first significant digits. 

Figure [4] shows results of numerical simulations (discrete 
points) versus theoretical predictions (solid curves) for a 5- 
state LZC-model. In this case, all coupling constants are as- 
sumed to have the same value, and probabilities were plotted 



as functions of this size of all couplings. Results appear to 
be in excellent agreement with the theoretical prediction. Fig- 
ures |4j a > c ) clearly show that the probability to remain in the 
initial level saturates at a nonzero value, which is independent 
of the coupling constants, as far as couplings are sufficiently 
large and all but 0-th level slopes have the same sign. 



III. COUNTERINTUITIVE BEHAVIOR IN 
NON-INTEGRABLE MODELS 

In this section, we explore the probability to stay at the ini- 
tially occupied level with non-Coulomb potentials. We con- 
sidered a set of diabatic potentials of the 0-th level of the form 



e (t) = q/t r , 



(22) 



with some constant parameters q and r. In particular, for 
our tests, we chose 5 -state models with q = 0.2 and 
r G {2, 3/2, 1/2}. All other parameters were chosen as in 
Fig. [4j a )> i.e. for the case when all other states have posi- 
tive slopes of the diabatic energies. Our results for transition 
probabilities are summarized in Fig. [5] 

Figure [5] shows that the phenomenon of saturation of the 
probability Po^o at large couplings is not an artifact of the 
integrability of the model. For example, Fig |5ja) shows that 
for the decay of the 0-th diabatic energy as 1 / y/t in a five- 
state model, the survival probability reaches a minimum at 
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FIG. 4. Numerical check of the transition probability formulas in 
(4+l)-LZC model for the case of equal coupling constants: gj = g, 
and slopes given by f3j = b 2 — b\ * (j — 2) for j £ {1, 2, 3, 4}, and 
(a) k 2 = 0.2, bi = 0.25, b 2 = 1.25; (b) k 2 = 0.2, bi = -0.2, 
b 2 = -0.65; (c) k 2 = 0.15, ft = -0.5, f3 2 = -1. Solid curves 
are theoretical predictions. Discrete points represent the results of the 
numerical calculations. Insets on the right illustrate time-dependence 
of diabatic energies. In all numerical tests, the time evolves from 
t = 0.00001 to t = 800 with time step dt = 0.00001. 



some coupling g between states, and then it is actually grow- 
ing. Figure |5[b) shows that the decay of 0-th level as l/£ 3 / 2 
produces qualitatively similar behavior to the integrable case 
at r = 1. For r = 2, the probability to stay at the 0-th level, 
in principle, slowly decays with increasing g but this decay 
is not comparable to what one would expect from the inde- 
pendent crossing approximation. Semiclassically, one would 
expect that the transition probability should be close to the one 
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FIG. 5. Numerically obtained transition probabilities in a (4+l)-state 
model with power-law decay of the diabatic energy of the 0-th level: 
e (t) = q/t r , for q = 0.2 (a) r=l/2, (b) r=3/2, (c) r=2. All coupling 
constants are equal: gj = g, and slopes are given by f3j = b 2 — b± * 
(j - 2) for j e {1, 2, 3, 4}, and bi = 0.25, b 2 = 1.25. Red curve 
is the semiclassical prediction of Eq. {23} based on the independent 
crossing approximation. 



predicted by the independent crossing approximation: 



0^0 



=1 _ 



(23) 



where bk 



-rfii is the slope of the 0-th level near the in- 



tersection of its diabatic energy with the diabatic energy of 
the i-th level. The red curve in Fig. [5] shows the prediction of 
the approximation ( [23] ). This approximation works relatively 
well up to g — 0.4 but then it hopelessly fails. For example, 



for g = 3.5 and r = 2, Eq. (23 ) predicts Pq 



10 



-43 



i.e. 



even for such a moderate coupling strength, the difference be- 



tween the prediction of Eq. ( 23 ) and the numerically obtained 
result is over 40 orders of magnitude. 
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IV. LZC HAMILTONIAN IN THE MODEL OF COUPLED 
QUBITS 



One may wonder if the LZC-model with more than two 
states can be realized experimentally. Here we show that a 
3-state LZC-model is actually found in the fundamental sys- 
tem of two coupled qubits. The latter is the minimal hardware 
block needed to build a quantum computer, which explains 
why it was realized in a number of solid state structures and 
in atomic gases. Our exactly solvable model can serve as a 
test for the accuracy of a control over such device structures. 

Our model is described by the Hamiltonian: 



H = J(t) (1 



a 1 



Z 2 ) 



-B z (t) (a] 



(24) 

where a l a is the a-component Pauli operator acting in the 
space of the i-th qubit. First term in (24) describes the 
Heisenberg-like coupling between qubits. Second and third 
terms describe magnetic-field-like couplings. These parame- 
ters are usually required to be controllable, i.e. time dependent 
on demand, for quantum information processing. 



the following choice of parameters is assumed: 



y/2' 



(25) 



with constants k and g. Here we note that one does not have to 
reproduce singularity at t = because at sufficiently strongly 
coupling J, transitions among states will be suppressed, so it 
is sufficient to keep the qubit system at the ground state with 
large J for t < and then continuously decouple them. In this 
case, the state ^ ( | + | |t) ) decouples from other three 
states: 



|o> = on) - 1 m .. 
|i) = ltt), 



|2> 



U)- 



(26) 



(27) 



(28) 



In the basis (|26|i-(|28|i the Hamiltonian ( |24| i has the following 
matrix form: 



H 



-k 2 /t g 
9 Pt 
9 



(29) 



Hamiltonian (24) would correspond to the LZC model if which corresponds to a 3-state LZC-model. 



V. TWO-STATE LZC MODEL 



A. Relation to confluent hypergeometric equation 

The two-state LZC model, 

d k 2 d 
i—a = — a + gb, i — b = f3rb + ga, (30) 
dr r dr 

is trivial to solve. Here we discuss some of the possible its applications, which to some degree, can be extended to multi-state 
LZC. 

By changing variables a — » ra' and t = r 2 e~ %7r / 2 f3, then differentiating 2nd equation in (30) over time, we find the confluent 
hypergeometric equation: 

d 2 

t-^b+(r-t)b + sb = 0, (31) 

where 

1 - k2 1 -fg 2 k 2 \ 

are constants. The solution of (31) that satisfies initial conditions b(0) = 1, a(0) = is given by the Kramer's function 
M(s, r, r). This solution has asymptotics 

M(*,r,0) = 1, M(s,r,r -+ oo) - w , (33) 

r» 

where f(r) is a time-dependent phase factor that does not influence the transition probability. Using that + ix)\ 2 = 
7r/cosh(7rx), we obtain the probability to stay in the initial level: 

p-^g 2 1 'ft _l p-^k 2 

/3>0: P ^o= Pi^i= - , , (34) 
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FIG. 6. Left: Probability to stay at level at time t —t +oo for f$ = 1 > 0. Solid curves are theoretical predictions by Eq. ( [34) . Dashed 
curve is the Landau-Zener prediction obtained by linearizing time-dependence of diabatic energies near the avoided crossing point. Discrete 
dots represent the results of the numerical calculations. Right: Probability to stay at level at time t — >> +oo for f3 = — 1 < 0. Solid curves 
are theoretical predictions by Eq. ([35]>. Discrete dots represent the results of the numerical calculations. 



< : P ^o = Pi-n = 



-^k 2 +g 2 /\0\) 



1 + e 



(35) 



B. Alternative to LZ-formula 



Equation ( |34| ) takes into account the nonlinearity of the crossing curves. Hence it can be used as a better approximation than 
the LZ-formula. Consider, for example, an arbitrary model of a two-state avoided crossing: 



id = e(t)a + gb, e(t) 
ib = f3tb + ga, 



-I3t + -t 2 



(36) 



where n is half of the radius of the curvature (in the energy space) of the diabatic energy near the curve crossing point. It is easy 
to verify that LZC model corresponds to k = k/(f3 3 / 2 ). Equation (34) shows that if k —> oo, i.e. if the level crossing is linear, 



LZC solution transfers into the standard Landau-Zener formula. Solution of LZC shows that when k is finite, the limit of large 
g and the adiabatic limit do not generally coincide. 

Figure [SjLe/O shows the effect of the curvature radius of the level. It clearly demonstrates that the Landau-Zener formula 
has only a limited range of validity, restricted to large k. The deviations from the Landau-Zener prediction become particularly 
strong at k 2 < 1. In fact, the simplicity of Eq. ( 34 ) and the fact that the transition probability in the LZC-model depends on two 



rather than one independent parameters suggest that Eq. ( |34| ) can be a better approximation at intermediate sizes of the curvature 
and the level coupling than the Landau-Zener formula. Figure ^Left) confirms our expectation for the case of quadratically 
decaying diabatic energy e = 1/t 2 . Although Eq. (34) is not much more complex than the Landau-Zener formula, Fig. [^Z^/iO 
shows that it represents a substantial improvement at intermediate coupling strengths. Figure ^Right) also shows results of the 
numerical check of Eq. ([35) at f3 < 0. 



C. A model of avoided diabatic level crossing 

If diagonal elements of the Hamiltonian do not intersect but pass close to each other, one can speak about nonadiabatic 
transitions at avoided diabatic level crossing. Two-state models of avoided diabatic level crossings usually can be expressed in 
the form 

id = e(t)a + gb, e(t) = eo H — t 2 + . . . , eo > 0, 

ib = ga. (37) 

The minimal set of important parameters for such models are the minimal distance between diabatic levels eo, the relative 
curvature k at the moment of reaching this distance (which characterizes the rate of divergence of levels), and the coupling 
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FIG. 7. Left: Comparison of approximations based on LZ and LZC models for a model of an avoided level crossing. Discrete points represent 
results of numerical calculations of the probability to stay at the level at time t — V +oo in the model of quadratic decay of 0-th diabatic 
level: eo(£) — q/t 2 , and f3 — 1/2 in Eq. (36 1. Solid and dashed curves are the theoretical predictions of, respectively, Eq. {34} for LZC 
and the Landau-Zener formula in which parameters are tuned to produce the same level crossing rate and curvature as defined in Eq. (36) . 
Right: Comparison of approximations based on LZC model and numerical results for the model with quadratic diabatic energy of the 0-th 
level: e(£) = eo + t 2 / (2k). For the latter model, discrete points represent results of numerical calculation of the probability to stay at level at 
time t — ► +oo if the evolution starts at t — ► — oo (not at t — 0). Solid line is the theoretical prediction of Eq. ( [35) for LZC which parameters 
are tuned to produce the same values of eo and k in Eq. (37). 




between the diabatic states g. In case of LZC, a simple phase transformation maps this model at f3 < to Eq. (37) with 
e(t) = \ fi\t + k 2 /t, as illustrated in Fig. [s| 

Parameters eo, k are straightforward to derive for the LZC model: 



2k 



(38) 



Eq. ( 38 ) can be inverted, 



*o/(2 



(39) 



Unfortunately, it seems that LZC model with parameters ( [39] ) cannot be used as a good approximation to an arbitrary model 
of avoided diabatic transitions with the same minimal distance and curvature because of a strongly skewed shape of the curve 
j3t + k 2 /t. Figure [TJ^/g/zO compares probabilities to remain at the 0-th level in a model with a quadratically changing diabatic 
energy and the LZC model with the same values of basic parameters provided by ( [39) . Reasonably good agreement is found 
only at relatively small values of k and g. 

Finally, we note that the theory of nonadiabatic transitions in two-state systems has been very well developed. LZC is only 
one of many useful solvable systems in this field CHS. Moreover, numerical simulations of such systems is a trivial exercise. 
The two-state LZC is a solvable model with a constant off-diagonal coupling, which is, sometimes, a required assumption. The 
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simplicity of its transition probability formulas makes the LZC model a useful alternative to the LZ-formula in cases when 
considerable accuracy is not required but rather a quick nonperturbative estimate of a transition probability is needed. 

I 



